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We revisit some of the recent neutrino observations and anomalies in the context of sterile neutri- 
nos. Among our aims is to understand more clearly some of the analytic implications of the current 
global neutrino fits from short baseline experiments. Of particular interest to us are the neutrino 
disappearance measurements from MINOS and the recent indications of a possibly non-vanishing 
angle, 6113, from T2K, MINOS and Double-CHOOZ. Based on a general parametrization motivated 
in the presence of sterile neutrinos, the consistency of the MINOS disappearance data with addi- 
tional sterile neutrinos is discussed. We also explore the implications of sterile neutrinos for the 
measurement of |C//i3| in this case. We then turn our attention to the study of |C/e3| extraction in 
electron neutrino disappearance and appearance measurements. In particular, we study the effects 
of some of the additional CP phases that appear when there are sterile neutrinos. We observe that 
the existence of sterile neutrinos may induce a significant modification of the &13 angle in neutrino 
appearance experiments like T2K and MINOS, over and above the ambiguities and degeneracies 
that are already present in 3-neutrino parameter extractions. There are reactor experiments, for 
instance those measuring disappearance like Double-CHOOZ, Daya Bay and RENO, where this 
modification is less significant and therefore the extracted \Us3\ value when sterile neutrinos are 
present is close to the one that would be obtained in the 3-neutrino case. Based on our study, we 
also conclude that the results from T2K imply a 90% C.L. lower-bound on \Ue3\, in the "3-1- 2" 
neutrino case, which is still within the sensitivity of future reactor neutrino experiments like Daya 
Bay, and consistent with the one-cr range of sin^ 2^13 recently reported by the Double-CHOOZ ex- 
periment. Finally, we argue that for the recently determined best-fit parameters, the results in the 
"3-1- 1" scenario would be very close to the medium/long baseline results obtained in the "3-1- 2" 
case analyzed in this work. 
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I. INTRODUCTION 

Neutrinos have now been unequivocally established to be massive particles, but with very small masses. Experiments 
over the past two decades have firmly established a framework of neutrino oscillations that describe solar, atmospheric 
and reactor neutrino experiments (see for instance [T] and references therein). 

The experiments are consistent with the existence of three electroweak eigenstates (i/g, z^^, 1/,-) and three mass- 
eigenstates (z^i, 1^2^ ^'s)- While the absolute neutrino mass scale has been very difficult to measure, the mass squared 
differences between the mass-eigenstates {Am^i, Am§2) are known to good accuracy [T}{3]. Also, two of the mixing 
angles in the lepton sector (6*12, 6*23) are known to good significance [2H1| and the third (^13) is being measured and 
will be measured to better and better accuracy by current and forthcoming experiments [SHS]. It is already clear, for 
instance, that the mixing in the lepton sector is very distinct from the quark sector. 

Inspite of these spectacular successes there are still many outstanding questions related to neutrinos. For instance 
it is not understood why neutrinos have such tiny masses or why their mixing angles are so much different from the 
quark sector. There have also been discrepancies from various short-baseline experiments that have been very hard 
to accommodate in the three active-neutrino picture. This has led to many studies incorporating additional singlet 
neutrino states to the framework 'Q'-TBl. 

Our main focus in this paper will be to gain a better analytical understanding of scenarios with additional singlet 
neutrino states and how they may affect current and forthcoming medium/long-baseline neutrino experiments. We 
are particularly interested in the measurements from MINOS and the determination of the reactor angle ^13 at various 
medium/long-baseline neutrino experiments. A study similar in spirit to ours was done in [191 , for the case of an 
additional sterile neutrino. The focus of our study will albeit be different from theirs and will also be motivated by 
the current global fits, incorporating two additional sterile neutrinos. 

In sections II and III we briefly review the current state of neutrino experiments and the viability of sterile neutrinos 
in the context of the standard model and cosmology. In section IV we briefly outline the short-baseline limit of neutrino 
oscillations and the global fits based on them. We also fix our notations here. In section V we study the implications of 
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these short-baseline global fits to medium/long-baseline neutrino experiments and explore various theoretical features. 
In section VI we give a summary of our results. 



II. CURRENT STATE OF NEUTRINO OBSERVATIONS 



Let us briefly review the current state of neutrino parameters. In a three neutrino framework let us label the v 
mass eigenstates by latin indices, % e (1, 2, 3), and electroweak eigenstates by greek indices, a e (e, /i, t). 

In the three neutrino framework, various experiments have measured the two mass-squared differences to be [IH3] 

[Amgal ~ 2.4 x 10"^ el/^ , (1) 
lAmiil ~ 7.6 X 10"''^ el/2 _ 

The overall mass-scale is not determined from oscillation experiments alone, but cosmological considerations imply |20l 

EI] 

^m^. < 0.6 &V . (2) 

i 

This still leaves an ambiguity in the ordering of the mass eigenstates. If Awgj > neutrinos are said to be in a 
normal mass hierarchy (NH) and if Arn|2 < they are said to be in an inverted mass hierarchy (IH). 

Similar to quarks and the CKM matrix, the electroweak and mass eigenstates in the lepton sector are related by 
a mixing matrix. The relevant angles in this Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix [52] are 
denoted by 6*23 (related to atmospheric oscillations), 613 (reactor oscillations), and 612 (relevant to solar oscillations). 
They are currently measured to be in the intervals [iHS] 



37° < 923 < 56° , 
0° < 6*13 < 13° , 

30° < 9i2 < 38° . (3) 

There are very strong indications [3] that the atmospheric angle 623 is close to maximal (7r/4). Solar data [3] also 
favors a 'Large Mixing Angle' (LMA) solution that indicates large 6*12. These intervals roughly translate to a PMNS 
matrix of the form 



7 /cxp. 
^PMNS 



0.9 0.5-0.6 

0.6 0.3- 0.7 0.6 - 0.8 | . (4) 
0.5 0.5 - 0.8 0.6 



Some comments are in order. The above matrix is very close to a tribimaximal mixing matrix which has vanishing ^13 
and maximal 9i2, 023- The question of CP-violation in the lepton sector, even in the simplest case of neutrinos being 
Dirac particles, is at present open. This is intimately tied to the question of a vanishing ^13, since in the event of a 
vanishing or extremely small value the relevant Jarlskog-invariant (Jcp area of the unitarity triangle oc sin 2^13) is 
also zero. 

Though almost all the observations from atmospheric, reactor and solar neutrino experiments can be accommodated 
conveniently in a three active-neutrino framework, there are a few outliers that may suggest existence of additional 
sterile neutrinos (i.e. neutrinos that are electroweak singlets). Let us briefly review them. 

LSND employed a liquid scintillation detector to search for De oscillations [53] with a baseline of about 30 m. 

The neutrino energies were such that L/E ^ 0{1) (Km/GeV). An excess of about 88 events was observed in the 
20 MeV ^ E < 60 MeV range. A conventional oscillation interpretation of the data requires Am^ ~ 0.1 - 10 eV^ and 
is referred to as the LSND anomaly [23 . 

MiniBooNE was designed to test the LSND anomaly. They searched for both z/^ — )■ and — De oscillations [24] . 
The E 1 GeV) and L (541m) were both much bigger than LSND, but L/E was still kept to be about the same. 
No oscillations in neutrino mode was observed above 475 MeV, but an unexplained 3cr excess of fe events was found 
below 475 MeV. Subsequently, a 2.8cr excess of i>e was observed in the range 475 MeV < E < 1250 MeV consistent 
with LSND. There are preliminary updates from the MiniBooNE collaboration for their anti- neutrino data [25] where 
there is an excess below 475 MeV similar to ly data, but the —?' Vf. signal above 475 MeV has now diminished from 
before [25]. 
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Recently, a re-evaluation of the reactor anti-neutrino flux [55] also indicates a deficit from the expected value. The 
data sets used for the evaluation came from various experiments with very short baselines (L < 100 m). The observed 
to expected flux ratio was found to be 0.943 ± 0.023 compared to a previous value of 0.979 ± 0.029. 

There was also a puzzling discrepancy between muon neutrino and anti-neutrino disappearance data from the 
MINOS experiment [21 [57] where the best-fit for the v data exhibited a higher Am^ and lower sin^ 2Q than the v data. 
Since it is a survival probability being measured, a resolution of both v and v data requires an effective CPT violation, 
say for instance through some non-standard matter interactions iBD]. There is an updated measurement from the 
MINOS collaboration in the v mode giving Am^ < 3.37 x lO'^ eV^ at 90%C.L. assuming sin^ 2^ = 1.0 12S]. The 
most current, preliminary best-fits give |Ato^| = [2.621°;^^ (stat) ± 0.09(syst.)] x lO'^eV^ sin^ 26* = 0.95l5J:i'j'(stat.)± 
0.01(syst.) for the new anti-neutrino data [28]. The previous discrepancy between v and v data therefore seems to 
be almost resolved. Prior to these measurements by MINOS, the strongest constraints on v parameters were from a 
global fit dominated by Super-Kamiokande data that included both atmospheric and i/^ results [551 US] • 

There have also been indications for a non vanishing ^13 from various experiments |31H33j . A non-zero dn has 
far-reaching implications for CP-violation in the lepton sector and the consistency of other neutrino parameter fits. 

MINOS detected Vf. appearance at the far detector [HP, compared to expected background, suggesting for Scv = 

2sin2 6l23sin2 26li3 < 0.12 (NH) , 

2 sin^ 023 sin^ 26*13 < 0.20 (IH) . (5) 

The best-fit points [M] for 2 sin^ 6*33 sin^ 26*13 are deduced to be 0.04li°:E;;^J for normal (NH) and 0.079i5J;|^53 for 
inverted hierarchies (IH). lAmj^l = (2.32;g J2) x lO'^ eV^ |Am^i| = (7.591°;^^) x lO'^ eV^ 6*33 = 0.785 ±0.001 and 
di2 = 0.60 ± 0.02 have been assumed in the above fits. 

The T2K experiment [35] observed six ly^ events that pass all selection criteria at the far detector. This suggests, 
at 90% C.L., again for Scp = 

0.03 < sin2 26li3< 0.28 (NH) , 

0.04 < sin^ 26li3< 0.34 (IH) . (6) 

The best-fit points [32 for sin2(26'i3) are found to be O.llljj j^g (NH) and 0.14^°:^^ (IH). The above T2K limits and 
best-fit values are extracted for sin^ 26*33 = 1-0 and Am§3 = 2.4 x 10^"^ eV^. 

Most recently, there is a preliminary result 33 from Double-CHOOZ, based on the first 100 days of data, showing 
at 68% C.L. 

sin^ 26*13 = 0.085 ± 0.029 (stat.) ± 0.042 (syst.) . (7) 

This result is particularly interesting since, being a reactor neutrino disappearance measurement, it is independent of 
CP phases and the mass hierarchy. 

Some of the above short-baseline discrepancies may be interpreted as being due to the presence of extra electro- weak 
singlet neutrinos. With this in mind, in the next section we briefiy discuss the viability of sterile neutrinos. 



III. STERILE NEUTRINOS 



From the invisible Z-decay width and LEP measurements, the number of active-neutrino species is constrained 
to [J 

^z-w.dth ^ 2.92 ± 0.05 , 

TV^'^^ = 2.984 ± 0.008 . (8) 

Also, as already mentioned, most of the current neutrino oscillation data can be accommodated in a three active- 
neutrino framework. The Troitzk [33] and Mainz [5S] experiments give a mass bound of about 

< 2.3 eV , (9) 

by measuring the endpoint region of the tritium /3-decay spectrum. In cosmology, neutrinos play a significant role by 
effecting the expansion history and growth of primordial perturbations which lead to a tighter mass bound [50] [5T] 

< 0.6 eV , (10) 
for three fiavor mixing. The KATRIN experiment (351 is speculated to reach a sensitivity of < 0.2 eV. 
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All these nevertheless still leave open the possibility of other neutrino species that are singlets under the SM 
gauge groups and therefore "sterile". Trying to resolve the LSND and MiniBooNE anomaly with data from solar 
and atmospheric neutrino measurements require, as we noted in the last section, Ato^^^^.j^ ^ 0(1) eV^. A similar 
mass squared difference is also required to reconcile the reactor anti- neutrino flux deficit. There are also some very 
controversial indications from the Heidelberg-Moscow experiment |37j of detecting a neutrino mass 0.17 eV < mpp. 

An analysis [38], a few years ago, combining data from cosmic microwave background (CMB), large scale structure 
(LSS) and Lyman-a constrained the mass of a fourth sterile neutrino to be rUg < 0.23 eV assuming they are thermal. 
In the non-thermal case it was shown that the constraints are non-trivial in the mass-density plane, but still viable. 
This was re-emphasized in a study [31] that combined the WMAP 7-year data, BEN, small-scale CMB observations 
and measurement of the Hubble parameter from the Hubble space telescope. The study concluded that the current 
data set mildly favors extra radiation in the universe and derived constraints on the number and mass of possible sterile 
neutrino species. More recently, in [?D], the authors conclude that though sterile neutrinos are disfavored by hot dark 
matter limits in minimal ACDM, extending the standard cosmological framework to include additional relativistic 
degrees of freedom or a dark-energy equation of state parameter w < —\ can relax these constraints substantially. 

Constraints were also put on active-sterile mixing by measuring neutral-current interactions at MINOS |41j . By 
measuring the depletion of the neutral current event rate at the far detector, a 90% limit was placed on the fraction 
of active neutrinos that transition to a sterile neutrino (assuming 0i3=O) |41j 

/. = . < 0.22 (0.40) . (11) 

The number in the bracket is for the assumption ^13 = 11.5° and 5cp = tt. 

Apart from experiments reviewed in the last section, we should also mention two other experiments - SAGE [42] 
and GALLEX [33] ■ They give for i/^ disappearance a measured to calculated ratio R — 0.86 ± 0.05 consistent with 
each other. If interpreted as due to an additional sterile state, this corresponds to a Am^j^^jj^ = 2.24 eV^ and 
sm^2dee — 0.50 [33] ■ This conclusion has been weakened though by a recent analysis [3S], based on its consistency 
with KARMEN and LSND data. 

Motivated by all these indications, extensive studies have been performed on fitting the short-baseline neutrino 
discrepancies to various models with sterile neutrinos and checking their consistency with other experiments [5141 3|. 

The inclusion of additional sterile neutrinos to the three active ones adds more structure to the neutrino oscillation 
formalism, without changing the basic framework. "3 -I- s" refers to the case of 3 active and s sterile neutrinos. The 
case of 3 active neutrinos in the standard model shall henceforth be denoted as 3z^SM. For 'n' neutrinos (active and 
sterile) the total number of angles in the mixing matrix is n(n — l)/2. The angles that rotate sterile states to sterile 
states are not relevant for neutrino oscillations and hence this number can be trimmed to 3(n — 2) angles. Similarly, 
the number of CP phases (Dirac) that could be present in the PMNS matrix is (n — 2)(7t, — l)/2. The number of 
physical CP phases relevant to electroweak physics and oscillations is 2n — 5, after field re-definitions. For "3 + 1" 
this gives 6 angles, 3 phases and for "3-1- 2" this gives 9 angles, 5 phases. 



IV. NEUTRINO OSCILLATIONS AND THE SHORT-BASELINE LIMIT. 

The probability for a neutrino state a oscillating into a state j3 in vacuum is given by 

Pii^a ^ >^p) = Sap - AJ2m*a^Up^UaJU*^^) sin' {Aml^ ^) -f 2 ^ 5([/^ C/^, f/^^- ) sin(Am2^.^) , (12) 

i>j i>j 



where i,j denote the mass eigenstates and a,l3 include both active and sterile neutrino eigenstates. 3? and 5 stand 

= Ami J 



for the real and imaginary parts respectively. We will also use the notation Aij = Amf^L/AE. 



In most oscillation experiments (which are constructed with sensitivity to a particular Am?) one can simplify the 
above by taking a two neutrino limit. In this limit we get the familiar result 



1 - sin^ 26>sin^(1.27 Am^(eV^) H^^^l ) ; a = /3 , 
sin' 20 sin^(1.27 Am^eY^) ^^^^^ ) ] a P ■ 

We note a few well known properties. In the two neutrino limit P^'^{va — > vp^Ueff) — P^^{vi3 — >■ Va',Ueff)- In 
general P{va — ^ i^p]U) — P{vp — >■ Va]U*). Under the assumption of CPT, in addition to having the mass-squared 
differences same for both neutrinos and anti-neutrinos, we also have P{va — > Vj3\U) = P{vp — Va'iU). The above 



two results lead to the fact that P{i'a — >■ i^p',U) = P(i'q — )■ vp]U*). Thus, note that it is the last term in Eq. (12) 
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Model 


Am|i(eV") \U,4\ 




Amii(eV2) |f/e5i 






" 3 + 2" 


0.47 0.128 


0.165 


0.87 0.138 


0.148 


1.64 



TABLE I. Global best-fit points using SBL data for the "3 + 2" case [TT]- 



that distinguishes neutrinos and anti-neutrinos (when the number of families is greater than 2) for a 7^ /3, indicating 
CP violation (CPV). 

It is also important to emphasize that for disappearance measurements or survival probability (a — /3) the last 
term vanishes. Hence, survival probabilities in vacuum will not be effected by any CP phases and must be the same 
for both neutrinos and anti-neutrinos if CPT holds. The last term being an odd-function of /S.m? is also sensitive to 
the mass hierarchy of the neutrino species in principle. 

We will be interested in two specific limits of Eq. ( [l2| . The first limit is the short-baseline (SBL) limit which is 
relevant approximately when L ~ 0(1) Km, E ^ 0(1) GeV and consequently L/E ^ 0(1) (Km/GeV). The other 
limit we would be interested in is the long-basehne (LBL) limit where L > O(IO^) Km and L/E > O(IO^) (Km/GeV). 
An intermediate case where L ~ 0(1) Km, E ^ O(10~^) GeV and L/E O(IO^) (Km/GeV) is often called medium- 
basehne (MBL). 

We are going to primarily focus on a " 3 -f 2" scenario as analyzed in the recent global fits [TT] , with additional mass 
squared differences in the 0(l)eV^ range. The comprehensive global fits, to SBL neutrino experiments, in '11] use 
appearance data from LSND MiniBooNE [24], KARMEN J46!, NOMAD ^47j along with disappearance data from 
Bugey [3S], CHOOZ [5], Palo Verde [15] and CDHS [SD]. The most recent analysis further includes full spectral data 
from SBL reactor experiments ROVNO [5T], Krasnoyarsk [55], ILL [S3] and Gosgen through rate measurements 
as summarized in [2S]. The main conclusion in the study is that the global fits to SBL oscillations, for a "'3-1- 2" 
case, improves significantly with the inclusion of the new reactor anti-neutrino flux data, though some tension remains 
in the overall fit [TT] . The global- fit values of [TI] from SBL experiments is shown in Table |l] We will adopt these 
values to analyze MBL/LBL measurements. Rather than taking the numbers in Table |l] as numbers set-in-stone, our 
attitude rather will be to view them as quantifying, to good extent, potential effects of sterile neutrinos in SBL. 

The SBL limit is the case most suitable for analyzing the LSND and MiniBoone experiments since they both were 
designed with a characteristic L/E ^ 1 (Km)/(GeV). A short-baseline (SBL) assumption leads to a simplification of 



the most general oscillation formula in Eq. (12 1. The following approximations may be made 



Am^2 -> 



Ami, ^ , (14) 

compared to Am'j, where i > 3. This is partly motivated by requirements from LSND and MiniBoone observations 
which require 

^'m.l^.na -0{l)cY'\ If there were no sterile neutrinos with AtoJ^^^^,^ ^ 0(1) eV^, then under the SBL 
approximation P{i'a ^ i^a', U) — 1, as should be expected since the baseline is not sufficient for significant oscillations 
into other flavors. 



Using the above assumptions in Eq. ( 12 ) we get for the "3-1- 2" case 



P„/3(CPV)sBL = 5a/3 " 43?(xL^)sin2A,i -43?(^„0)sin2A54+ ^ 23(xL^) sin 2A,i + 23(Ca;3) sin 2A54 , (15) 

i=4,5 i=4,5 

where 

i^fi = {U*r^Up5Uc.4U;^) . (16) 

For a — j3 — e this gives, 

P(Pe ^ i>e)sB.= 1 - 4(1 - |C/e4p - | f/eS H [peil^ siu^ A41 + |f/e5p siu^ A51] - 4 1 f/^g H C/e4 | ' siu^ A54 . (17) 

When Q ^ /3, the CP phase(s) can play a major role and the appearance or conversion probabilities may be different 
for neutrinos and anti- neutrinos. It may be shown that in the SBL approximation for "3-1-2" there is only one relevant 
phase [TT]. Consider the case a — ^ and /3 = e. Define the relevant phase as (also see Appendix A) 



6 = aTg{u:,u^4Ue5u;,) 



(18) 
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Neutrino/Anti-Neutrino survival probabilities (SBL) 




Neutrino Energy (GeV) 



Conversion probabilities (SBL, 6=1 .64jt) 




Energy {GeV) 



FIG. 1. Survival {Pee ) and conversion (P^e) probabilities for neutrinos and anti-neutrinos in a SBL experiment with L = 541 m, 
assuming S = 1.647r and parameters from Table [I] The assumed distance to the detector corresponds to the baseline in 
MiniBooNE |24| . The difference in neutrino and anti- neutrino conversion probabilities is clearly visible in the bottom figure. 



Then 



mf.e) = \U:^U^,Ue5U;,\cOsS, 



(19) 



^ixte) = -\Ue,\'\U,^\' - \U:,U,,U,5U;,\ COS 6, -sixlJ = \U:,U^^U,5U;,\ sin S , 

^ixle) = -\Ue5\^\U^5\' - \U:,U^iUe5U;,\ COS S, 3(x^J = - | C/;4 t/e5 f/;5 I ^ • 

As an aside, it should be mentioned that in the corresponding " 3 + 1" case there is no CP phase that is relevant to 
the SBL limit and therefore the conversion probabilities in that case cannot distinguish between v and D, for instance, 
at MiniBooNE. 



With the inclusion of a CP phase, from Eq. (15), we have for the conversion probabilities 

P{i^ii i^e;CPV)sBL = P{i^t, Ve]5 = 0)sBL + 8|^^e| sin A41 sin A51 cos(A54 -5)- COSA54 



(20) 



P(f^ i/e;CPV)sBL = ^'(t'M Ve]S = 0)sBL + 8|^^e|sinA4isinA5i cos(A54 + 5) ~ cosA 



154 



(21) 
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To give a more holistic picture we plot the disappearance and appearance probabilities in Fig.[T] assuming Mini- 
BooNE base length, with the matrix elements and CPV phase of Table |l] We see that for low energies the depletion 
of is greater than v^. in the survival probability. Also note the expected enhancements and attenuations in the 
conversion probabilities, in this energy range, for D and v. The presence of a non zero CP phase could therefore, in 
principle, explain the difference between v and D measurements in MiniBooNE |55| . 

As we shall see, the effects of sterile states at the near detector (ND) will have implications also in LBL measurements 
since usually it is a ratio between the far detector (FD) and ND neutrino fluxes that are compared, to probe for 
oscillations. Thus a depletion of the neutrino flux at the ND due to sterile states will affect MBL/LBL measurements 
too. 



V. LONG-BASELINE LIMIT OF NEUTRINO OSCILLATIONS. 



We now turn our attention to LBL experiments and it is our intention to understand in more detail the consequences 
of the SBL fits, or equivalently the presence of the two sterile neutrino states, to these experiments. Specifically, we 
explore in this section the impact of the SBL global fits to MINOS disappearance measurements and the recent ^13 
determinations. The questions we would like to address are whether MINOS disappearance data can constraint or 
discriminate in a definite way the " 3 + 2" case from 3i^SM and if the recent indications of a significant, non- vanishing 
reactor angle ^13 are effected drastically by sterile neutrino states, if they exist. 

In the LBL limit the assumption is that L/E (Km/GeV) may be anywhere in the range 10'^ — 10* or higher. The 
larger values are more appropriate for atmospheric and solar oscillations. For very large Am^, the oscillations get 
averaged out inside the detector and we may replace sin^(Am^2^) by its expectation value (1/2). Consequently, 
Am|]^ and Am|]^ are now averaged out since they are 0(1) eV^. Am2i may also be set to zero to first approximation. 
We will include this to be non-zero later, when discussing ^31 determination. 

Calculating the oscillation probability from Eq. ([T2|, in the above limit, for "3 + 2" gives 



Pa/3(CPV)^BL = 5o.p ~ 23i{xc,p) - m^c^p) sin^ A54 - 43fi(Ca/3) sin^ A32 + 5(Xa/3) + 23(ea/j) sin 2A54 

+ 23(Ca/j)sin2A32 , (22) 

where 

Xa/3 = {5al3 " U^iUp^ - UasUpr-,) ^ U*iUpi = ^ xL/3 : 

i=4,5 i=4,5 

C,al3 = (SaP - UasUps - UaiU^i - t^aS ^^^5) (^a3 '^/JS ) • (23) 

For a survival probability {a ~ j3) this may be re-written as 



^ (1 - 25ft[xo 



sm A32 



1 - 25ft[Xaa] 



(24) 



Here we have neglected the term proportional to ICq/sI which tends to be negligible numerically, since it is quartic in 
the sterile neutrino matrix elements. 



We observe from Eq. ( 24 ) that in the case of additional sterile states there is an effective normalization factor 
modifying the survival probability as well as a modification of the coefficient (^ sin^ 29 in 3z^SM) of the energy 
dependent term. Thus, we conclude that the LBL survival probability in the presence of sterile neutrinos is of 

a general form 

Pl^l^ ^ AA„ [1 - sin^ 2,9f ■ sin2 A32] , (25) 

with 

sin^ 2df- - . (26) 

1 — 2v 

The normalization factor (A/'q), in principle, may be extracted by looking at asymptotically large neutrino energies. 
The sin^ 2z?°^' coefficient determines the dip of the first minima. The modification to sin^ 2'd, when there are sterile 
neutrinos, has two parts - a direct modification of Caa and a further scaling of this quantity by Afa- Observe that the 
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quantities appearing in the LBL limit, Afa and the modifying terms in sin^ 2^^' , are completely determined by SBL 
measurements through the matrix elements. 

Using the global fit values from Table |lj the various quantities that appear in the survival probability, Eq. (22), may 
be computed as (for an assumed Ue3 — 0.1 and [/^s = 0.707) 



Xee 


= 0.0342 


(3i/SM 


0), 




= 0.0467 


(3i/SM 


0), 


Cee 


= 0.0095 


(3i/SM 


0.0099) 




= 0.225 


(3J/SM 


: 0.25) . 



The numbers in the brackets are the corresponding values in the 3i^SM case. Using the above values we get 

p3+2 ^ 0.932 [1 - 0.0408 sin^ A32] , 
= 0.907 [1 - 0.993 sin^ A32] , 



(27) 



(28) 



which is relevant to any FD in an MBL/LBL experiment, for instance, the MINOS FD. For comparison, in the 3i/SM 
the expressions corresponding to above would have been 



P: 



ee 



= 1 - 0.0396 sin^A32 , 
2 



- i - sm A32 



(29) 



Compared to the 3:/SM prediction the overall flux is reduced by about 7% for and by as large as 10% for 
in Eq. (28). The effective angle is seen to be not modified significantly in the survival probability and may still, in 
principle, be extracted to yield a value that is close to the true value. We will come back to this point again while 



discussing |?7^3| and \Ue'i\ determination. From Eqs. (27) and (28) we note explicitly that the modification to sin^ 2i9, 
when there are sterile neutrinos, is due to a direct modification of Qaa and a further scaling of this quantity by Ma- 
lt is also to be re-emphasized that when the sterile neutrino AtoJ^^^;,^ is large there will also generally be an effect 
in the ND, apart from the above effects in the FD, especially for low E. Since neutrino experiments usually compare 
a ratio of the fiuxes at the FD and ND, taking into account geometric and other effects, modifications at the ND due 
to sterile states may also become relevant along with the FD effects. We will in fact see that for the values in Table 
|T] the ND effects at MINOS are not completely negligible. 

Although we specifically derived the expressions and numerical values for the " 3 + 2" case, a similar parametrization 
should he valid for any "3 + s" scenario with AtoJ^^^;,^ ~ 0(l)eV^. Note also that owing to Eq. (25) being a survival 



probability, any CP phases that may be present are completely irrelevant. This is particularly important since the 
CP phases that are relevant in LBL experiments may in general be mutually exclusive to the ones that are relevant in 
SBL experiments 11]. Thus, the above features must be applicable independently of CP phase structures in the LBL 
limit. The CP phases relevant to LBL may nevertheless become important in neutrino appearance measurements as 
we shall see. 



A. Sterile neutrinos and the MINOS disappearance data. 

We now turn our attention to the analysis of MINOS muon neutrino and anti- neutrino disappearance data | 
It is interesting to ask if the MINOS disappearance data can discriminate or put constraints on scenarios with sterile 
neutrinos, specifically "3 + 2", as motivated by the SBL global fits. The quantity of interest here is the ratio of the 
observed events at the MINOS far detector (L" = 734 Km, 4.2 Kt fiducial mass) to that expected at the far detector 
if there were no neutrino oscillations. The latter is extrapolated from the MINOS near detector (L"™°^ = 1.04 Km, 
23.7 1 fiducial mass) through a Monte-Carlo, taking into account geometric and pion kinematic effects [51 [?fl 1^ . 

Simple fits to the MINOS neutrino (7.25 x 10^" protons on target (POT)) and anti-neutrino (2.95 x 10^° POT) 
data [51 1271 [25] are shown in Fig. [2] with various parametrizations, motivated by Eq. (25). For the anti-neutrino analysis 
we have taken the new preliminary data i28j . 

The values for the fit parameters are shown in Table [ll] along with their la errors. In the first parametrization (Case 
I) both Ato^ and sin^ are floated, as in the analysis of the MINOS collaboration. A physical boundary on sin^ 2'd 
is not imposed in the fit. In the second parametrization (Case II) Am^ and the normalization are free parameters 
with sin^ 2i?^ fixed at 1. The implicit motivation here is that, as we saw previously, the sin^ 2i?p coefficient is modified 
only minimally in disappearance measurements even when sterile neutrinos are present. In the final parametrization 
(Case III) we fioat A/'p, sin^ 2^^ and Ato^. In this context, it is worth re-emphasizing that the MINOS data shown 



is actually the ratio of the FD and ND neutrino fiuxes and there will be deviations from Eq. (25), especially at low 
energies for the values given in Table |l) due to ND effects. 
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Muon Neutrino events at MINOS far detector. 



1-sm'^(2())sin^(A ra'^ L/4E): sin^(29)=0.953,A ra =0.0023 (x?n.d.f: 1.84) 



N(1-sin'(28)sin''(A L/4E)): 



, N(1-sin^(A nf L/4E)): N=0.892, Am,=0.00206 (xjn.d.f: 
^ L/4E)): N=0.895,sin^(2ft)=0.946 ,Am^=0.0021 6 (yt/n.d.f: 



Theoretical FpQ/Fp^Q (x /n.d.f: 



1.67) 
1.67 
1.53) 



O 
6 



CO 
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Muon Anti-Neutrino events at MINOS far detector. 



1-sin^(2())sin^(Ajn^ L/4E): sin^(2O)=C.86,Ajn^=0.0C3 (x^/n.d.f: 1.24) 

N(1-sin^(A L/4E)): N=1.16, Am,=0.00305 (x,/n.d.f: 1.61) 

N(1-sin^(2d)sin^(A L/4E)): N=1 .07,sin^(2ft)=0.883 ,Am =0.00305 (xp/n.d.f: 1.33) 

Theoretical Ppq/F^q (x /n.d.f: 3.73) 




5 10 15 20 

Anti-Neutrino Energy (GeV) 



FIG. 2. Fits to MINOS neutrino 3 and anti-neutrino [271 128] far-detector data assuming three different parametrizations 
motivated by Eq. ( 25 I - Case I (red) , Case II (dark blue) and Case III (purple) . The theoretical curve from Eq. ( 30 1 (light blue) , 
calculated based on the SBL gobal-fit values, is also shown. Though the plots are shown only till 20 GeV, the fits were done for 
the full range of neutrino energies in the data set (50 GeV for both and v^). The neutrino data corresponds to 7.25 x 10^" 
POT and the anti-neutrino data to 2.95 x 10^° POT. 
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Case 




sin^ 2i9^ 




xVn.d.f 


I 


2.31 ±0.10 


0.953 ±0.04 


It 


1.65 


II 


2.07 ±0.09 


It 


0.895 ±0.03 


1.48 


III 


2.17 ±0.13 


0.946 ± 0.048 


0.897 ±0.03 


1.48 










1.53 



Case Am2(10"^eV^ 



sin^ 2i9, 



xVn.d.f 



I 

II 

III 



3.0 ±0.23 
3.05 ±0.2 
3.05 ±0.22 



0.86 ±0.08 it 1.24 

it 1.16 ±0.13 1.61 

0.883 ± 0.086 1.07 ±0.122 1.33 

3.73 



TABLE II. Parametric fits to MINOS neutrino and ant-neutrino data 3, 27, 28 for various parametrizations that take into 



account the possibility of sterile neutrinos, motivated by Eq. (25 1. The symbol f denotes that the particular parameter is not 
floated in the fit under consideration. 



The theoretical ratio of the flux at the FD and ND may be estimated from Eqs. (15 1 and (22) as 

Fpo(a*) ^, 1 - 2x^p - 4C^p sin2(AAf|2%^ 



MINOS 



'^MM 17 l' \ — n „ r MINOS ■ ["^^ ) 

This ratio estimated from the SBL global-fit values is also plotted in Fig. [2] |Am§2| — 2.32 x 10^^ eV^ has been 
assumed to calculate the theoretical curve. In the MINOS data set, the ND flux has been extrapolated to the FD, 
assuming no oscillations, including effects from beam-line geometry and meson decay kinematics [51 1271 [^ . i?^^ and 
Rfifi, which are obtained assuming naively a coUimated neutrino beam, may be directly compared to the real data 
sets to good approximation due to this. To quantify how the ratio predictions compare to the data sets, we may 
calculate a reduced-x^ (x^/n.d.f). It is found that the theoretical ratio prediction has a reduced-x^ of 1.53 and 3.73 
for the neutrino and anti-neutrino data sets respectively (Table |ll| . 

From the neutrino data we extract an overall normalization, 0.897 ± 0.03, that is close to the one theoretically 



calculated in Eq. (28) from the SBL global fits, which had a central value of 0.907. The extracted sin^ (0.946 ± 
0.048) is slightly lower than that predicted from the SBL global- fit values (0.993), but still within la. The uncertainties 
are larger in the anti-neutrino data set and the number of data points is also smaller. For this case it is found that 
both the extracted normalization and sin^ 2i? are significantly larger and smaller respectively than that predicted 
from SBL fits. There is hence some potential tension of the global-fit values with anti-neutrino data. Note that since 
this is a survival probability, if CPT is satisfied the probabilities in vacuum should be the same for neutrinos and 
anti-neutrinos. 

In Fig. [3] we plot the 68% and 90% likelihood contours for our simple fitting procedure, as applied to the MINOS 
neutrino data. Note that systematic uncertainties have not been taken into account as nuisance fit parameters. With 
the inclusion of a possible normalization J\f^ in the fits, the extracted Am^ decresases. The 90%C.L. lower bound on 
sin^ 2-3 also shifts marginally to a smaller value (0.90 — ?► 0.886) as we include Af^. Although we fix A/"^ in determining 
the 68% and 90% contours in Case III, the conclusions don't change significantly even if we include an error bar on 

In the plots and fits shown, we have not changed the bin size from that given in the publicly available data set 
(the fits as done by the MINOS collaboration are for a bin size that is different from ours). We have nevertheless 
checked that re-binning the data (consistent with Fig. 2 of [3]) does not significantly change our observations. The 
MINOS collaboration also perform a more comprehensive likelihood analysis, compared to our fitting procedure, 
incorporating systematic uncertainties as nuisance parameters. For comparison, the values of the best-fits from 
the MINOS collaboration are \Am'^\ = 2.32l°;J| x lO'^ eV^ sin^ 2z9 > 0.90 (90% C.L.) for neutrinos |3] and \Am^\ = 

[2.62l°;^^(stat) ± 0.09(syst.)] x lO-^eV^ sin^ 2d = 0.95;[! }?(stat.)±0.01(syst.) for the new anti-neutrino data [HIISH]. 
So, although the central value is slightly different from ours, the lower bound on sin^ 2i? is very close to the one obtained 
in our fit. The comparison between the 3i/SM and "3 + 2" cases shows only a small shift of this lower bound and 
we expect this to be the case even after systematic errors are included. Thus our simplified analysis and observations 
are still meaningful as long as we are comparing parameters and confidence levels extracted using the same fitting 
procedure, as in Table |ll] and Fig.[3j 

Let us now try to derive some simple scaling relations between the 3i^SM and " 3 + 2" confidence intervals. In 
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Likelihood Contours (MINOS, Case I) 
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FIG. 3. Likelihood contours from our fitting procedure for MINOS neutrino data. The two plots are for Case I (left) and Case 
III (right) as defined in Table |ll] Systematic errors have not been included. It is observed that our estimate of the 90%C.L. 
bound for sin^ 2i9 shifts from 0.885 to 0.853, along with a reduction in Am^, when we include a normalization factor Af^. 



general, from Eq. (24), it may be shown that if a is a lower-bound (say the 90%C.L. limit) on the coefRcient of 
sin^ A32 in P^^^, the matrix element |C/^3| virould be constrained as 



\U, 



m3 \app. 



(31) 



when the 'incorrect' assumption of 3;/SM and no sterile neutrinos is made. With the assumption as " 3 + 2" we would 
have a constraint on the "true" value 



\U, 



IJ.3\"tr 



where 



/3; = l-|C/^4p-|f/^5| 



= 1 - 2x^^ = 1 - 2^; (1 - 131) 



(32) 



(33) 
(34) 



It is easy to sec from the above equations that the allowed interval for the extracted \ 1 is shifted slightly, to lower 
values, when one includes sterile neutrinos in the extraction procedure. Also note that in the case of the sterile- 
neutrinos, due to the modification, the coefficient o/sin^ A32 (which would have been ^ sin^ 2^23 *^ the ZvSM case) 
can never be unity. 

For the sin^ 2-d > 0.885 (90% C.L.) Si^SM (Case I) bound of Fig.|3] (left plot) this translates to 

It^A^al^pp. e [0.33,0.67] , (35) 

and for the or the sin^ 2?? > 0.853 (90% C.L.) "3 + 2" scenario (Case III), Fig.js] (right plot), it becomes 

IC^msIw e [0.29,0.66] . (36) 

The extracted |?7^3| indeed shifts downwards by a few percent from the 3z/SM case, when "3 + 2" is assumed. This 
will be relevant to us when we discuss Ve appearance measurements, since the value assumed for |J7p3| will have a 
bearing on the extracted [C/esj in that case. 

It is clear from the analysis of the MINOS neutrino and anti-neutrino data sets that they by themselves cannot 
distinguish between the 3j^SM and "3-1- 2" cases in a definite way. There nevertheless is a marginal improvement 



12 



in the reduced-x^ , in the neutrino fits, when sterile neutrinos are included. The normalization as extracted from 
the neutrino data is close to that theoretically predicted for the FD and including any ND effects from sterile states 
improves the fit slightly. The uncertainties in the anti-neutrino data are larger and the reduced-x^ obtained for the 
theoretical prediction is poor. With higher statistics the analysis in the anti-neutrino sector may be improved in the 
future. 

We will now explore the effects of sterile neutrinos on 613 determination, in the context of disappearance and 
appearance measurements. 



B. Effects of sterile neutrinos on the determination of ^13 



There have been recent indications from experiments like T2K, MINOS and Double-CHOOZ for a possibly non-zero 
013 . Let us try to understand the implications of sterile neutrinos for ^13 determination in these present and upcoming 
MBL/LBL experiments. If there are sterile neutrinos with AmJ^^^jj^ ^ C'(l) eV^, then the mixing angles corresponding 
to them may be extracted in principle solely from SBL experiments. These matrix elements thus extracted, may then 
be used in the analysis of MBL/LBL experiments that aim to measure the reactor angle ^13. Due to the presence 
of the additional sterile states there may be possible modifications to the extracted 813 angle or more precisely the 
extracted matrix element | Ues \ ■ 



1. Reactor i/^ disappearance measurements 



Let us first consider neutrino/anti-neutrino disappearance experiments which measure survival probabilities. The 
CHOOZ [S] experiment and its successors Double-CHOOZ "S^, Daya Bay [7] and RENO fH| are examples in this 
category. We take Double-CHOOZ as an example for purposes of our discussion. At Double-CHOOZ the FD is at 
1050 m and the ND is expected to be placed 400 m from the reactor cores. For these base lengths and E 3MeV) 
we may use the MBL approximation at the FD. 

We have for the Ve/Vf, survival probability, at the FD 

- 1 - 2 (|[/,4P + |C/e5p - P.it - If/esi' " Pe^U,^?) - 4 (|C/,3p(l " \U,:i?) - (1 - /^DlC^eSp) siu^ A31 

- 4 (|C/e2|'(l - |C/e2n " (1 " P'e)\Ue2?) siu^ A21 + 8 1 C/e2 H (JeS I ' siu A21 sin A31 COS A32 , (37) 

where, as before 

/?^ = l-|[/e4p-|;7e5r (38) 

Note that we have now retained Am^i terms explicitly. Since this is a survival probability, any CP phase that 
may be present in the MBL/LBL limit is completely irrelevant and does not cause any ambiguities. This is not the 
case, as we shall see, when we consider conversion probabilities where the phases may play a significant role. For 
the Double-CHOOZ baseline and characteristic E 3MeV), the terms which are proportional to A21 and quartic 
in sterile-neutrino matrix elements may be dropped without incurring significant errors. This leads to the familiar 
expression 



derived m Eq. (24) before. This may now be rewritten as 



where 



~l-2Xee-4CeeSin2 A32, 

sin2 2^, sin^ A32] , 



(39) 
(40) 

(41) 

(42) 
(43) 



AAe = 1 - 2Xee ^ I - 2{l ~ /?^) , 

sin^2.. = ^^.4MtK). 

1 - 2Xee 1 - 2(1 - /30 

Under the 3-neutrino assumption this survival probability may have been written as 



^ 1 - 4|C/e3|^,p.(l - Pe3\Lp) A 



32 



From this we may associate an "apparent" value for the ^13 angle, through the relation 
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l^eslapp. 4|f/e3|;pp.(l-|C/e3|'pp.) \Ue3\ 4 | j7e3 P (1 - | ^e3 1" ) 



0.224 0.19 0.219 0.183 

0.198 0.15 0.194 0.145 

0.168 0.11 0.165 0.106 

0.133 0.07 0.131 0.067 

0.087 0.03 0.085 0.029 



TABLE III. Comparison of some representative apparent and true values of measured | Ues \ and 4 | Ues P (1 — | U^s in a neutrino 
disappearance experiment. A " 3 + 2" scenario is assumed and the matrix elements for the estimate are taken from the SBL 
global fits. 



leading to the standard form for the Ve ~^ i^e survival probability 

P-,^l-sin2 2^^^f^■sin2A32. (44) 
In general, if a' is an upper-bound on the coefficient of sin^ A32 in the survival probability P^e, we have the constraint 



2 ^ 1 - \/l - a' 



\Ue3\Up. ^ ' (45) 

for the apparent value when the "incorrect" assumption of no sterile neutrinos is made. With the correct assumption 
we would have a constraint on the true value 



f3' - J/3'l - Ka' 

|C7,3|Le ^ \UeZ? < . (46) 



Compared to the {U^^l extraction case in Eq. (32), we have dropped a solution that is already ruled out. Also note 
that in the probability expressions, leading to me above extraction, the A/'e factors cancel while taking the FD/ND 
ratio. This has to do with the fact that both the FD and ND probability expressions are usually more akin to the 
MBL limit in reactor neutrino experiments, for instance for a ND around 400 m as in Double-CHOOZ. 



In Table III we list for comparison, values of |J7e3| and 4 |C/e3p(l — |?7e3p) determined in the neutrino disappearance 



case, using the 3i^SM ("apparent") and "3 + 2" ("true") assumptions. The true values for jj/esl are found to be 



generally smaller than the apparent values, as is to be expected from Eq. (46 1. Nevertheless we observe that, due to 
the absence of CP phase ambiguities, negligible matter-effects and smallness of the sterile-neutrino matrix elements, 
the value of lU^sl extracted in reactor neutrino disappearance experiments are not modified significantly even when 
sterile neutrinos are present. Due to this, the extracted value of ^13 even under the incorrect assumptions is still close 
to the true value. 

The above observations are especially pertinent in view of the recent preliminary result from Double-CHOOZ [33], 
suggesting 

sin^ 29^3'=''°°^ = 0.085 ± 0.029 (stat.) ± 0.042 (syst.) . (47) 
With a Double-CHOOZ FD at 1050 m and a future ND at 400 m, the FD/ND ratio will be of the form - MBL/MBL 



and as pointed out before Eq. ( 46 ) the normalization factor Afe would cancel in the numerator and denominator in 
this case. Although Double-CHOOZ currently lacks a ND, they normalize their measurement to the Bugey experi- 
mental data [H], which accounts for an approximate 0.945 suppression factor with respect to the expected one. The 
result is then fitted to the form 1 — sin^ 29i3 sin^ A32 to extract the value of sin^ 2^13 . The difference between the 
normalization factor extracted from the Bugey experiment and the Afg predicted in the "3-1- 2" neutrino scenario is 
small compared to the systematic and statistical errors. We expect that this small difference, between the energy 
independent normalization factors, will not lead to any relevant variation of the extracted sin^ 2'd value obtained from 
the fit. Of course if the full energy dependence for events at the Bugey detector at 15 m is considered, one should use 
an expression TZ^e analogous to 7^^;^ used in Eq. (30) for MINOS. 



Taking the central value above and assuming that there are two additional sterile states, characterized by the values 
in Table |lj gives 

\Ue3\ir°"^ - 0.147, 4|C/e3|L.(l - | C/e3 lL. = 0.085 , 

|f^^^|D-cHooz ^ ^ 4|i7e3|'(l - |C/e3n°-°''°°" = 0.082 . (48) 

As emphasized before, the difference between the extracted values in the 3-neutrino scenario and in the "3-1- 2" 
scenario is very small compared to the current errors. 
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Now, consider a standard parametrization of the " 3 + 2" PMNS matrix 

3 

^'i.L= n (49) 

j>i,i—l 

where the product is to be done from right to left and Mjj is a complex or real rotation matrix in the ij-plane (see 
Appendix A). Let us generically denote by 9^ the small sterile angles {Oij, j > 3). Using this parametrization, it may 



be seen that the coefficient in Eq. (41) is still very close to sin^ 26*13, the deviations being of 0{6g) (Appendix A). 
The quantity |J7e3| is nothing but sin^ia in the 3i^SM case. In the "3 + 2" case, using the standard parametrization, 
it gets modified to cos 6*14 cos ^15 sin ^13 (see Appendix A). The quantity 4|C/e3p(l — |J7e3p) in the 3i^SM case would 
have corresponded exactly to sin 2^13. In the "3 + 2" case it deviates from sin 26*13 by terms of 0(6*2) (Appendix 
A). Note also that Jcp oc sin 26*13 in the 3i'SM case. 

All the above conclusions are also applicable to upcoming experiments like Daya Bay [7J and RENO 8 that aim 
to measure ^13 through a disappearance measurement. This is to be contrasted with appearance measurements that 
we discuss next. 



2. Ue appearance measurements 



The value of ^13 may also be deduced by looking for v^. appearance in experiments that measure neutrino conversion 
probabilities. This includes experiments such as T2K 132] and MINOS |5T] . 

Let us focus on T2K [32] for the purposes of our discussion. At T2K the neutrino energy peaks around 0.6 GeV 
and the ND(s) and FD are at 280 m and 295 Km respectively. This ensures A32 ~ tt/2 at the FD when E ^ 0.6 GeV, 
giving an oscillation maximum. These base lengths also enable us to use all the approximations for LBL again at the 
FD. The ND effects can again be quantified and studied as before. 

The transition probability in vacuum, keeping Am^i explicitly, may be written as 

PIb^ - 4|C/^3|'|C/e3|' sin^ A31 + 4|C/^2nC/e2p sin^ A21 + 8|C/;3||?7e3||[/^2||C/:2| sin A31 sin A21 cos(A32 - 83) 

+ 4|C/,,3l|f/e3ll^1 sin A31 sin(A3i - <5i) + 4|C/^2||C/e2||/3l sin A21 sin(A2i - ^2) + 2(|C/^4p|(7e4p + If/^sl'l^ 

+ |;7^4||C/:4l|(:^;5l|C^e5|cOS<5) , (50) 



where 



|3" = J2K^Ue^. 

i>4 



S, = arg (C/;3C/e3/3") , 
<52 = arg {u;,U,2f^") , 

<53 = arg {u;3Ue3U,2U:,) . (51) 



The only approximation we have made in Eq. ( 50 ) is to average terms containing large sterile mass-squared differ- 
ences. In contrast to the previous case there now appears non-trivial CP phases, as this is a conversion probability. 
Due to the presence of these phases there could be interesting interferences between the various terms and it's seen 
that one can no longer drop terms, as we did in the survival probability case, without significant errors. 

These phases which appear in the LBL limit are in general independent of the phase S extracted from the SBL 
global fits. Also note that when /3" the phase is not independent and is given by 6i — 62- The effective CP 
phases Si and S2 may be related to the 'fundamental' CP phases ^12 and 613 in some particular parametrization of 
^3+2^ (Appendix A). 



To get a better understanding of what the various terms in Eq. (50 1 mean, we briefiy look at the corresponding 



expression in the 3i^SM case. In the 3:^SM the conversion probability has the well-known form 



^ ^ItT + i^l"'" + 2 JP|-"P3,-r C0S(A32 - <53) , (52) 



with 



P^^f^' ^ sin^ 923 sin^ 2^13 sin^ A31 , 

Pl'^^" ^ cos^ 6*23 sin^ 2012 sin^ A21 . (53) 
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The first and second terms in Eq. (52 1 are the atmospheric and solar oscillation contributions. The last term denotes 
an 'interference' between the atmospheric and solar oscillations with a relative phase shift S3. Label this term Pq^Tatm 
and in terms of the matrix elements it is 

P™_"^f,f. EE 8|C/;3||f/e3||C/^2||f/e*2l sin A31 sin A21 cos(A32 - S3) . (54) 

If we define Sep — — arg([/e3), then for small |J7e3| values, ^3 in the Si^SM is almost equal to —Sep- The P^'^^" term is 
0{A2i) and small for most experiments we are interested in. The P^J^^tm term is superficially sensitive to the mass 
hierarchy, since under +|Am32| — )■ — lAm^jl it picks up a negative sign and the argument (IATO32I— (S3) — >■ (|Am32|+(S3). 
Nevertheless, it should be noted that there is no actual sensitivity in the vacuum case, since a rescaling of the CP 
phase, S3 TT — S3, would undo the above transformation 56J. 

For later comparison to "3 + 2", in Fig.|4]we make some illustrative plots in the 3i^SM case, assuming T2K baseline 
(295 Km). Note that in the 3z/SM case 

8 1 (7^3 n C/e3 P = 2 sin2 023 sin^ 2f?i3 . (55) 

Plotting this combined quantity in Fig.|4] (top left) and later allows us to readily consider a non-maximal atmospheric 
sector in \Ue3\ extraction. The CP phase S3 now is almost equal to —Sep, for small |t7e3|. Note from Fig. [i] (top right) 
that, for a fixed value of the CP phase, the theoretical differences between NH and IH can be more pronounced as 
we move away from A32 ~ 7r/2. This of course does not imply any actual sensitivity to the mass hierarchy in these 
measurements, due to the invariancc under S3 ^ tt ^ S3 and +|Am|2| — >■ ~|Am|2| mentioned earlier 1561 . 

The bi-probability plot, bottom figure in Fig.|4j shows the probability orbits in the (P^e, Pfis) plane. The orbits 
are traced as we vary S3, whose values may be read off from the color wheel at the origin, and the size of the ellipses 
are determined by the magnitude of \Ue3\- Since A32 Tr/2 the cosi53 contribution in the interference term is small 
and the ellipses get squeezed as A32 n/2, tending towards a line [ST]. Due to this there is no (^3,013) degeneracy 
in the plotted orbits. In this case, if matter matter effects are small (at T2K, for instance, they are small to good 
approximation), CP-violation can in principle be measured directly by comparing P^e and P^g- 



The first three terms in Eq. (50 1 corres pon d to these terms above - atmospheric, solar and an interference term 

ATM. 



between them. The next two terms in Eq. ( 50 1 are contributions solely from sterile neutrinos and modulate the Pf 



and P^^" with relative phase shifts of Si and S2 respectively. Let us label these terms APf^M.'^^ and APq 

APf;;7^i = 4|C/^3||C/e3||/3"|sinA3isin(A3i -(5i) , 

AP|"-^^ = 4|[/^2||C/e2||/3"|sinA2isin(A2i -52) ■ (56) 

The last three terms are energy independent residues obtained after averaging out terms involving large sterile 
neutrino mass-squared differences. Let us call these residue terms collectively as p^'^'^' , 

EE 2(|C/^4nC/e4p + |C/^5nC/e5|' + | C/^ 1 1 f/e*4 1 1 C^;5 1 1 ^^eS I COS <5) . (57) 

Though higher in order, they become important for small values of \Ue3\ to give a positive-definite conversion 

probability. It is also worth emphasizing that both Pq^Tatm ^'^d AP^ti^"^^ are oc |C/e3|- This observation will become 
relevant later when we try to understand cancellations among them. 

In Fig. [5] we show contour plots of the conversion probability P^e, for two fixed values of \Ue3\- T2K base length 
and characteristic neutrino energy have again been assumed. For the T2K ND distances (280 m) and E^, = 0.6 GeV, 
the effect of sterile neutrinos on the ND fiuxes is relatively minimal, but still leads to noticeable spectral distortion. 
We have included this effect in the analysis and the plots. Focusing primarily on ^ 0.6 GeV may be justified by 
the fact that the J-PARC v^^ beam has a very narrow side band [32] and in addition, under ideal conditions, most of 
the statistical power in sin^ 2(?i3 extraction may be expected to come from the region of the first oscillation maximum 
(A32 ~ 7''/2), tuned at 0.6 GeV. We will look at the effects of varying E later for comparison. 

Couple of things may be noted immediately from these contour plots in the {Si, plane, at A32 ^ 7r/2. The 
largest conversion to takes place in the vicinity of (0, — 7r/2) in both cases. In fact this is found to be true for all 

intervening values of \Ue3\ as well. In this region both P^^^tm and APfxiZ''^ have the same sign and constructively 
interfere with P^tm"- Heuristically, let us denote this situation as 

(0, -7r/2) : pA3.-V2 ^ p3.sM ^ TpiNT^^s^ ^ ^pP^^-^A , (58) 



where © denotes constructive interference and denotes destructive interference. As we decrease |?7e3| the conversion 
probability decreases as expected but the maximal conversion region is relatively unchanged. 
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FIG. 4. Plots in the St^SM case, assuming T2K baseline, showing extracted 8|C/^3p|f7e3|'^ for E — 0.6 GeV (top-left plot) and 
conversion probabilities in the energy bins {0.4GeV, 0.6GeV, 1.2 GeV} for fixed 8|C/^3|^|i7e3p (top-right plot). The thick-lines 
denote NH and the dotted lines denote IH. The top-left plot is for four fixed probabilities - 0.07, 0.05, 0.03 and 0.01. For 
comparison, the CHOOZ limit at 0.19 (90% C.L.) is shown. To avoid clutter, we have not shown in the plots the best-fit values 
from the other experiments - 0.04ll[];°^^ (MINOS-NH), Orndto^oll (MINOS-IH), 0.111E;;o6 (T2K-NH), 0.14l°;J^ (T2K-IH) 

vs. Pue) for five different values of 



and 0.085 ± 0.051 (Double-CHOOZ preliminary) [3TH33| . The bi-probability plot {Pp. 
4|[7e3p(l - |(7e3|^), correspoudiug to Table 



III 



is also plotted for the NH case. The color coding on them denotes the values of 
63 and the values may be ascertained from the color wheel displayed at the origin. Apart from \Ue3\, for simplicity, wherever 
applicable all other matrix elements have been assumed to be close to their tribimaximal values. 



Similarly, the lowest conversion probabilities occur in the vicinity of (ivr, — 7r/2), where -P™-Tatm ^^"^ APi^'-r^ °^ have 



the same sign but now destructively interfere with 



ATM. ; 



ATM. 



^ATM. ^ 



— ATM ^ ATM 



(59) 



Probably even more interesting is the observation of a thin band near 62 = 7r/2 for which the conversion probabihty 
bant ov< 
terms, 



is almost constant over the full range of i5i. In this band there is an almost perfect cancellation between the PI^atm 



and APfTM.''^ 



iySi,+n/2) : p^/2~-/2 



P. 



3l^SM , 



p. 



INT. — 63 



— ATM. 



(60) 
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FIG. 5. Contour plots of — >■ Ve conversion probability, with T2K parameters, as a function of two independent phases 5\ 
and 82 near A32 ~ ii/2 for the NH case. The T2K ND effects have been included. The color-coding denotes the magnitude of 
the conversion probability in each case. The quantity 4|[/^3|^(1 - It/^^sP) in the 3z/SM case would have corresponded exactly to 
sin^ 2^13. In the "3 + 2" case, as we had commented earlier, it deviates from sin^ 2^13 by terms of 0{9l). In terms of |L''e3| the 
above two plots correspond to 0.22 (left) and 0.085 (right), capturing the \Ue'i\ range in Table [TiT] In this case, besides \Ue-i\ 
and the global-fit values of Table |lj again all other matrix elements have been assumed to be close to their tribimaximal values. 



The conversion probability in this case should almost be identical to the conversion probabihty in the 3i^SM case, 

^3) ~ ±7r/2 in Eq. 
±7r. The conversion probabilities are 



when the interference term P^^ativ? in the 3z^SM limit almost completely vanishes (i.e. (A32 



(52)). In the 3i^SM case, near A32 ~ 7r/2, this happens at (Js = and ^3 



indeed found to match as expected upon comparison. Again, for the global-fit and close-to-tribimaximal values we 
are working with, this conclusion is seen to he true, independent of Uez , since both PqZ'atm '^'^^ IS.Patm. o,re cx | Ue3 \ ■ 



When 62 = 0, it is seen that the terms P. 



o_ATM. s-i^d APfrpM^''^ are out of phase with each other by 7r/2 in the phase 
Si (or now equivalently S3). Due to this, the conversion probability attains its maximum and minimum values at 
S3 — 7r/4 and ^3 = — 37r/4 respectively when A32 ^ n/2. 

In Fig.[6]we look at the probability orbits in the (P^e, Pp-s) plane for A32 ~ n/2. We see a rich behavior in the orbits 

depending on the value of S2- These may again be understood in terms of interference between P, 



— ATM 



3 ^p™'r-'53N 
\^ — ATM./ 



and APf;-'^ (APf: 



in P^e iPp.e)- Since P(i/„ iyp;U) 



P{va — J> vp]U*), all the phases change sign as we go 
from neutrinos to anti-neutrinos. Specifically as seen from Eq. (|60l, for S2 — —tt/2, the P0_atm 



and Apf^j.^ ^ terms 



cancel each other for P^g leaving a residue almost independent of 53. This may be symbolically expressed as 



^ — ATM. ^ ATM. 



(61) 



It is also interesting to note that for S2 ^ Q the (53,6'i3) degeneracy is reintroduced even when A32 ^ n/2. This 
is clearly seen from the fact that the elliptic orbits for adjacent |J7e3| values intersect. At the points of intersection 
both Ppg and P^e have the same magnitudes for different values of S3 and |J7e3|- This degeneracy can now mix 
CP-conserving and CP-violating solutions. An example of this may be observed in the ^2 = case of Fig. [6] where the 
orbit labelled by 0.183 intersects the orbit labelled by 0.106. The former solution is CP-violating {S3 ^ 0) whereas 
the latter is CP-conserving (S3 ~ 0). Note that in the bi-probability discussions, to first approximation, the effects of 
the base length (L) and neutrino energy {E) only appear through A32 ~ 7r/2. Due to this the general features of the 
probability orbits should be more widely valid, as long as we are in the vicinity of a conversion maximum. 

Fig. [7] shows the extracted values of 8|C/^3p|C/e3p for fixed conversion probabilities, again assuming that most of the 
statistical significance is coming from the region near A32 ~ 7r/2. The thick- lines indicate NH and the dotted-lines 
indicate IH for comparison. 
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FIG. 6. Probability orbits in the "3 + 2" neutrino scenario for the NH case. The respective values of 4|i7e3|^(l — It^esp) are 
labelled on the orbits. As before, apart from {Uesl ^-ud the global-fits, all other matrix elements have been assumed to be close 
to their tribimaximal values. It is observed that in contrast to the 3!^SM case, the (^3,^13) degeneracy may be re-introduced 
depending on the value of the phase 62- Also observe that ^2 ~ — 7r/2 and S2 = +n/2 are orthogonal choices for which Ppe 
and P^e remain almost constant respectively, as ^3 traces the orbit. These may be compared to the 3uSM probability orbits in 
Fig.|4] The value of 83 at any point in the orbit may again be deduced from the color wheel at the origin. 



The different values of 8|[/^3p|C/e3p extracted in Fig. [7] are easily understood by looking at the corresponding 
conversion probabilities near that particular CP phase region. For a fixed conversion probability, positive (negative) 

interference terms from Pq]:^^ and AP^t^j^^^ must be compensated by lower (higher) value of \Ue3\- Therefore, 
the variation of the extracted {Uesl niagnitude as we vary S3 must be anti-correlated with the conversion probability 
variation. 

Another point we would like to emphasize is that the actual \Ue3\ value extracted depends on the assumed value 
of \Ufj.3\, which may be extracted, as we saw in Eq. (32), from a disappearance measurement such as MINOS. We 



saw in the case of LBL, and specifically MINOS, that the presence of sterile neutrinos cause the value of extracted 
|J7p3| also to shift to lower values by a few percent 

\U,3\T. < |c/M3lit^ . (62) 
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FIG. 7. Extracted values of 8|C/;.3p|C/e3p in the "3 + 2" case, with T2K parameters {L = 295 Km, E = 0.6GeV), assuming 
fixed conversion probabilities - 0.07, 0.05, 0.03, and 0.01. The thick-lines are for NH and the dotted-lines are for IH. We 
re-emphasize that the quantity 8 1 (7^31^1(7631^ plotted would exactly correspond to 2sin'^ 6*23 sin^ 2013 in the 3vSM. The other 
matrix elements have been chosen as in the previous cases. For comparison, the preliminary best-fit value from Double-CHOOZ 
is at 0.085 ±0.051 |33|. 



Due to these considerations, plotting 8|(7^3p|(7e3p in Fig. fallows us to incorporate a non-maximal atmospheric sector 
and I (7^3 1 readily. Let us try to understand some of the theoretical features in the plots. From Fig.|4]we observe that 
in the Si^SM case the extracted |C7e3| is the same for normal and inverted hierarchy when 63 = ±7r/2V£'. This is 
due to the fact that the interference term becomes oc ±sin^ A31 for 63 — ±n/2. It is also clear that the maximum 



difference between normal and inverted hierarchies occur at — {0, ±7r}Vi?, as expected from Eq. (52 1. It is worth 
re-emphasizing that these theoretical features again do not imply actual experimental sensitivities to neutrino mass 
hierarchy, in vacuum, due to the invariance under 83 ^ t: — S^. 



Along similar lines we can understand the features in the "3-1- 2" case. The NH and IH cases should give the same 
extracted \Ue3\ when the terms which transform under jArngj] — >■ — jArngj], specifically P^^^tm ^^'^ ^Pa™.^'' , sum 
to the same numerical value, apart from the |(7e3| factor that is common to both. If a term by term equivalence is 
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FIG. 8. Conversion probabilities P^e in the " 3+2" case, for three different neutrino energies - 0.4 GeV (red), 0.6 GeV (green) and 
1.2 GeV (blue). The thick-lines are for NH and the dotted-lines are for Iff as before. They are plotted for a fixed 8|C/^3p|C/e3|^ 
of 0.085. It is clear that there may be significant differences between NH and IH depending on E and the CP phase structure. 



demanded (which is a stronger condition than required) between NH and IH, we must have 

[C0S(|A32| - (53) + C0S(|A32| + 5^)] ^ , 

[sin(|A3i| + 5i) - sin(|A3i| - <5i)] ^ , (63) 
These give the solutions in the (^i,(52) space : 

NH = IHV(L, i;) : (0, ±7r/2), (tt, ±7r/2), (-tt, ±7r/2) . (64) 
These regions are clearly visible in the top-left [82 = ^ 7r/2) and bottom {62 = +7r/2) plots of Fig.[7j with the 



understanding that (^3 = 5i — 62 

■~in 4- n 4- 4-1-11-^ n 1-1+ cm -t-ir-i /-\+ 1 

— ATM. 



Imposing the weaker condition that the net sum of ^q^Tatm ^^^^ APatm/^ be equivalent in the NH and IH cases 



leads to the solution 



NH . IH V (L, E) : . 2|^,2||^2l sin A^. cos IA32I 

^ ' cos((5i-52) |/3"|cos|A3i| ^ ' 
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This is independent of \Ue3\ and |C/^3|. For E = 0.6 GeV (A32 - n/2 for Lt2k), {Su52) ~ (7r/4,0) and (-37r/4,0) are 
among the approximate solutions to the above equation. This equivalence between NH and IH can be seen clearly in 
the top-right plot of Fig. [7] near these regions. 

For most of our theoretical discussions till this point, we were focused on the region near A32 ~ 7r/2 at T2K, and 
hence at neutrino energies near 600 MeV. As perviously noted, it may be argued that this is not too egregious a 
choice since the J-PARC beam has a very narrow side band and in addition, naively, most of the statistical power 
in sin^ 2^13 extraction must come from the region of the first oscillation maximum, tuned at 0.6 GeV. Nevertheless 
it is important to explore the variations with E, especially considering that the T2K appearance measurement 
observed 4 events outside the 0.6 GeV bin. 

Fig. Is] shows the variation in the conversion probability P^e as we vary E for NH and IH. They are plotted for a 
fixed 8lC/^3p|;7e3p of 0.085. We focus on the 400 MeV, 600 MeV and 1200 MeV energy bins. These bins correspond to 
A32 ~ 37r/4, A32 ^ 7r/2 and A32 ^ n/A for the T2K baseline. As is clear, the variations may be substantial between 
NH and IH as we move away from the oscillation maximum depending on {61,62)- We could now pose the question 
- for a fixed conversion probability near A32 ~ 7r/2 what is the smallest \Ue3\ it may be associated with for any CP 
phase structure in the Si^SM and " 3 + 2" cases ? 

In Fig.|9]we show the minimum attainable \Ue3\ values, in terms of 4|[/e3p(l — jC^sP), in the Sz^SM and "3 + 2" 
scenarios. The range of conversion probabilities in Fig. |9] correspond to those in Fig. [7] Note that the smallest possible 
\Ue3\ in the "3 + 2" case is always significantly smaller than that possible in the ivSM, 

\Ue3?ZlZ. > |t/e3|?lL. • (66) 

Also observe that, though reduced, the smallest \Ue3\ values in "3 + 2" are still different from zero, albeit extremely 

tiny for small conversion probabilities. This is because, Pq and APq by themselves are not sufficient to fulfill a 
given conversion probability in the range shown. 

In the 3-neutrino scenario, from the T2K lower-bound curves [32] (which give the 90% C.L. lower bound on sin^ 2^13 
as 0.03 — 0.04 for = 0), the smallest possible sin^ 2^13 is about 0.02 — 0.025, at 6cp — -~n/2. If we assume that 
the lower-bound curves correspond approximately to constant conversion probabilities, then using those approximate 
Pf_ie values in the "3 -f 2" scenario we may estimate, from Fig.joj a lower limit for 4|{/e3p(l — |i7e3p)- This gives for 
the " 3 + 2" case, 

4|C/e3p(l - iC/esP) > 0.008- 0.01; QUesl > 0.04- 0.05) at 90% C.L. . (67) 

Based on the above discussion, we note that in the "3 + 2" case, the results from T2K imply a 90% C.L. lower bound 
that is still within the reach (sensitivity) of future reactor neutrino experiments like Daya Bay [7, , and consistent with 
the la range of sin^ 26*13 recently reported by the Double-CHOOZ experiment. 



3. Matter effects 



All the above effects, due to the presence of sterile neutrinos, may be further modified by matter effects, depending 
on the base length and E. In the 3i^SM case these matter effects may be quantified as [58] 



^3.SM ^ ^ _ As V r'^^yr'' - sin2 2012 



2 sin (^M - l)A3i 

(A, - 1)2 
2 sin2(AM - l)A3i 



sin^ AmA3i 



''^'^ (^M - 1)2 

+ sin^ 26*1202 



2e si3 sin 29i2 sin 2^23 cos(A3i - 6cp) 



sinylMA3i sin(AM - l)A3i 



A^. 



A., - 1 



sin^ A^Agi 

23 A 2 



where 



_ Am^i 2EV, 

€ — — TT- , A 



Amii 



Eq. ( 68 1 is written to second order in e and sin ^13 , assuming a constant matter-density potential 

-14 ( Peru,, 



Kt ~ 7.56 X 10" 



Y. eV 



(68) 



(69) 



(70) 
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FIG. 9. Minimum possible \Ue3\ in the iuSM and "3 + 2" cases (NH and IH) expressed in terms of 4|l7e3|^(l - IC/es]^). In the 
plot, E is fixed at 0.6 GeV, the implicit assumption being that most of the statistical power in |(7e3| extraction may come from 
the vicinity of _E ~ 0.6 GeV (equivalently A32 ~ n/2 for T2K baseline) which is the oscillation maximum. It is clear that the 
lower-bounds on the extracted |(7e3| are generally much smaller with a "3 + 2" assumption, but still non-zero. The best-fit 
values in the ivSM for comparison are - 0.04lto osl (MINOS-NH), Qm^tl'lH (MINOS-IH), O.lltlJ og (T2K-NH), 0.14 
(T2K-IH) and 0.085 ± 0.051 (Double-CHOOZ preliminary) pTH33] . 



+0.11 

0.08 



where Pciust is the crust matter density and is the number of electrons per nucleon. For earth matter Ye — 0.5 to 
very good approximation. 

For a constant earth-crust density p^ruBt — 3 g/cm'^, we can estimate using the characteristic experimental parameters 
that 



^D-CHOOZ_ 0.0003, 

A™ ~ 0.06, 

^MINOS _ 0.3 . 



(71) 



Using the above values and Eq. (68 1 we can make estimates to convince ourselves that for Double-CHOOZ the matter 
effects are almost completely irrelevant and the extracted iJ/esI is hardly affected. For T2K, it is seen that the matter 
effects are still relatively minimal (near E ^ 0.6 GeV) but induce at most a few percent change in the extracted \Ue-i\ 
relative to the vacuum assumption. In MINOS the matter effects can become more significant and may induce larger 
modifications of the extracted \Uf.^\ somewhat obscuring any possible additional effects due to sterile neutrinos. 

In both cases above, for a fixed conversion probability, the effect of matter interactions is to decrease {increase) the 
extracted |J7e3| for NH (IH). This is probably most easily understood in a 2-neutrino limit by noting that the effect of 
the matter potential is to increase (decrease) the effective sin^ 26'm coefficient for NH (IH). Crudely, to lowest order, 
the above conclusion should still hold approximately in the "3-1-2" neutrino case. Also note that in contrast to matter 
effects, the direction in which |J7e3| was modified due to sterile neutrinos depended intricately on the (5i, ^2) CP phase 
structure. A comprehensive analysis of ^13 extraction including matter-effects at MBL/LBL, in the presence of two 
sterile neutrinos, is beyond the scope of the present work (see [33] and references therein in this context, for a " 3 + I" 
SBL fit incorporating matter effects). 



4- Comparison to the "3 -|- 1" case. 



It is interesting to point out in the MBL/LBL limit that for P^e: if one were to use the best-fit values in "3-1- 1" for 
the matrix elements, as in [12] say, the numerical values of /3" and the energy independent residue term p"^'^' defined 
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in Eq. (571 comes out to be numerically almost the same. In the "3 + 2" case we have 



1/3" 1-3+2" 



0.0351 




0.0021 



(72) 



and in the "3 + 1" case we get for the equivalent values, 




0.0354 



0.0025 



(73) 



These quantities along with the 2 independent phases are the only relevant quantities in the LBL/MBL limit that 
depend on the presence of sterile neutrinos. The larger Ar7i|j^ one obtains from SBL fits in "3 + 1", relative to "3 + 2" 
Am^ values, is irrelevant for MBL/LBL since the terms containing it get averaged at the FD anyway. 



In the present study, we revisited some of the recent neutrino observations in the context of sterile neutrinos and 
the global fits from SBL experiments, to understand their impact on current and upcoming MBL/LBL measurements. 

We noted that in general, for LBL experiments, the existence of sterile neutrinos lead to a distinct parametrization 
of the oscillation survival probabilities in terms of a normalization factor and a modified coefficient of the energy 
dependent term. We analyzed the MINOS neutrino and anti-neutrino disappearance data [3l |27l [28] from this 
perspective. Though the parametrization does lead to a marginal improvement in fit, it was found that the current 
MINOS data by itself does not definitively discriminate the " 3 + 2" scenario or the parameter values obtained from 
SBL fits. It was found that the jC/^al confidence interval shifts to lower values by a few percent when the possible 
existence of sterile neutrinos are taken into account. 

It was also commented that the recent measurements of a possibly non vanishing reactor angle may be affected 
by the existence of sterile neutrinos. We pointed out that the existence of sterile neutrinos may induce a modification 
of this angle (more precisely If/esl) in experiments that look at neutrino conversion probabilities, such as T2K and 
MINOS, and the perceived value may be shifted significantly from the "true" value in these cases. We also studied in 
detail the effects of additional sterile neutrino terms and their interference due to CP phases, in the "3 + 2" conversion 
probabilities. The probability orbits in the bi-probability plots also exhibited interesting features distinct from 3i^SM. 
It was, for instance, observed that the {63,613) degeneracy may be re-introduced depending on the CP phase structure 
in the " 3 + 2" scenario and that there may be orbits where either the P^g or P^e value remains almost constant with 
changing 63. 

It was also reiterated in the study that in the reactor experiments, these modifications due to sterile states are 
less significant. Due to this, the matrix element \Ue3\ when determined from survival probabilities under the 3:^SM 
assumption, is close to the " 3 + 2" value, as compared to when determined from conversion probabilities. Neutrino 
disappearance experiments include Double-CHOOZ |^ and upcoming experiments such as Daya Bay [7] and RENO [5] 
that will measure 6*13 to high precision. In this context we also conclude from our study that the results from T2K 
imply a 90% C.L. lower-bound on \Ue3\, in the "3 + 2" neutrino case, which is still within the sensitivity of future 
reactor neutrino experiments like Daya Bay [7] , and consistent with the one-cr range of sin^ 2^13 recently reported by 
the Double-CHOOZ experiment. Finally, we argued that the results in the "3 + 1" scenario, using the recent best-fit 
values, would be very close to the medium/long baseline results we obtained in the " 3 + 2" case. This was attributed 
to the numerical equivalence of the relevant parameters in both cases. 

Our analysis suggests that if the SBL global fits, including the anomalies, are in fact legitimate indications of 
sterile neutrinos in nature, then there may be interesting effects in MBL/LBL neutrinos experiments. A more 
comprehensive study in the "3 + 2" scenario including matter effects at MBL/LBL terrestrial neutrino experiments, 
such as MINOS [3T] and No^A [6^, is left for future exploration. We also plan to pursue in future, a study of how 
the "3 + 2" scenario affects survival probabilities for solar neutrinos [61j . using current global fit parameters. 

Note added : After the submission of this article a similar work appeared that discusses the recent Double- 
CHOOZ results considering the energy dependence of the events at the ND induced by the presence of sterile neutrinos. 
The final conclusion is similar to ours namely that the value of sin^ 2^13 is not significantly modified with respect to 
the three neutrino case, although due to the sterile neutrino effects the final uncertainties associated with the result 
are somewhat larger than the ones quoted by the Double-CHOOZ experiment. 



VI. SUMMARY AND CONCLUSIONS 
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Appendix A: Matrix elements 



For completeness, we list some of the relevant matrix elements in terms of the angles and phases in a standard 
parametrization, 

3 

The multiplication of matrices is to be perfomed from right to left. The rotation matrices may be real or complex. We 
choose a CP phase parametrization that is consistent with that employed in Under this convention the matrices 
Ri2, Ki3, 11^15, 1^34 and IR35 carry CP phases. 

The sterile neutrino matrix elements in this convention are 

9i4 , (A2) 

15 , (A3) 

025 sin 6*24 - e"^^= sin 614 sin ^15 sin 625 , (A4) 

925 . (A5) 

The sterile- neutrino angles (^14, ^15, ^24, ^25) and phase S15 can in principle be extracted from SBL measurements 
and specifically the global fits of Table |l] 

The active-neutrino matrix elements pertinent to our study come out to be 

Ue2 — cos 613 cos 9i4 COS ^156""^^^ sin 6*12 , (A6) 
Ue3 = cos 6*14 cos ei^e''^^' sin 6*13 , (A7) 

Up,2 — 008 6*12 cos 6*23 COS 024 COS 025 + C"*''^^ siu 012 (cOS ^13 (— COS 6*25 Siu 0i4 sin 024 — 008 0146"^^^ sin 0i5 sin 025) 

- cos 024 cos 0256"^'' siu 0i3 sin 023) , (A8) 

Up,3 = cos 013 cos 024 cos 025 sin 023 — e""'^^ sin 013 (cos 025 sin 014 sin 024 + COS 0146*'^^^ sin 015 sin 025) . (A9) 

The effective phases Si and 62 can be related in principle to the 'fundamental' CP phases 5i2 and 5i3 using the above 
relations. 





= cos 6 


15 


Ue5 




S 




— cos 6 
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f/M5 


— cos 6 
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